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Abstract. The high metallicity of the intra-cluster medium (ICM) is generally interpreted on the base of the 
galactic wind scenario for elliptical galaxies. In this framework, we develop a toy-model to follow the chemical 
evolution of the ICM, formulated in analogy to chemical models for individual galaxies. The model computes 
the galaxy formation history (GFH) of cluster galaxies, connecting the final luminosity function (LF) to the 
corresponding metal enrichment history of the ICM. The observed LF can be reproduced with a smooth, Madau- 
plot like GFH peaking at z ~ 1 — 2, plus a "burst" of formation of dwarf galaxies at high redshift. 
The model is used to test the response of the predicted metal content and abundance evolution of the ICM to 
varying input galactic models. The chemical enrichment is computed from "galactic yields" based on models of 
elliptical galaxies with a variable initial mass function (IMF), favouring the formation of massive stars at high 
redshift and/or in more massive galaxies. For a given final galactic luminosity, these model ellipticals eject into 
the ICM a larger quantity of gas and of metals than do standard models based on the Salpeter IMF. 
However, a scenario in which the IMF varies with redshift as a consequence of the effect of the the cosmic 
background temperature on the Jeans mass scale, appears to be too mild to account for the observed metal 
production in clusters. The high iron-mass-to-luminosity-ratio of the ICM can be reproduced only by assuming 
a more dramatic variation of the typical stellar mass, in line with other recent findings. The mass in the wind- 
ejected gas is predicted to exceed the mass in galaxies by a factor of 1.5-2 and to constitute roughly half of the 
intra-cluster gas. 
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1. Introduction 

The popular galactic wind (GW) sc enario for elliptical 
galaxies, introduced bv lLarsonl l|l974j) to account for their 
photometric properties, predicted as a side effect the pol- 
lution of the intra-cluster medium (ICM) with the chem- 
ical elements produced and exp elled by individual galax- 
ies llLarson fc Dinerstem|ll975|) . Met als in the hot ICM 
were in fact detected so on afterwards l|Mitchell et alll976t 
ISerlemitsos et al.lll977|) . 

Iron is generally used as tracer of the overall metal- 
licity, being the best measured element in the hot ICM. 
Typ ical iron abundances in the ICM are around 0.2-0.3 so- 
lar toenzinilll997tlFukazawa et al1ll998tlMushotzkvlll999t 
iFinoguenov et alJl2000t l200l|h The ICM seems to be also 
rich in a-elements, for which ASCA provid ed more firm 
estimates; data by iMushotzkv et al.1 l|l996l) yield [3A ~ 
0.2 ± 0.3 or, in the case of oxygen, [£] ~ 0.48±g;§^, 
which means, considering the uncertaint y, a marginal over- 
abundance of oxygen relative to iron. IMushotzkv et all 
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analysis bv iFukazawa et alT 1)1998'). who further suggest 
that [Jcj may increase with clus ter richness /temperat ure. 
Ilshimaru fe Arimotol l|l997l) and lArimoto et all l)l997|) ar- 
gue that considering the uncertainties and the revised 
value of the solar iron abundance, abundance ratios are 
consistent with solar. The first XMM and Chandra stud- 
ies seem to indicate the same l|Buotell2002|) . 

Recent studies have revealed a more complex distribu- 
tion of metals in the ICM. Gradients in t he iron abundance 
have been detected wi th ASCA data l|Finoguenov et all 
2000 j200lUWhitel2 000). a nd confirmed by BeppoSax dat a 
l|lrwin fc Bregmanl 120011 be Grandi fe Molendil l200lj) : 
in particular, sharp metallicity peaks in the cen- 
tral regi on seem to be typi c al of cooling flows /cD 
clusters jFukazawa et al l l200flt llrwin fc Bregmanl l200ll 
iDe Grandi fc Molendil l200l|) . [q/Fe] ratios also seem to 
display a composite behaviour, increasing with cluster ra- 
dius from ~ solar in the central regions to supersolar, 
typica l type II supernova (S N), values in the outer re- 
gions ijFinoguenov et alJ [2000). The existence of gradients 
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of metallicity and of abundance ratios is being confirmed 
by t he first XMM and Chandra results on a few clus- 



ters jlKaastra et all20_0lllTamura et allgOOll [Smith et alJ 



12003 llohnstone et all2002UEttori et alJ2002l) . buTno ex- 
tensive cl uster samples have yet been analyzed with these 
satellites l|Buotel l2002j) . It has been suggested recently 
that, in the presence of abundance gradients, emissivity- 
weighted estimates of the average metallicity might be 
higher tha n the true mass-weighted average, up to a fac- 
tor of two l|Pellegrini fc CiottfeOOlJ) . We shall not discuss 
abundance gradients in this paper however, but only the 
global average chemical evolution of the ICM. 

The source of such a large amount of metals in the 
ICM could be galactic, as in the original prediction by 
Larson, or reside in Population III pr e-galactic objects 
( White fc Reesll978tlLoewensteinl20"0ll) . The distinct cor- 
relation between the iron mass in the ICM and the lumi- 
nosity of elliptical and SO galaxies, 
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shown by lArnaud et alJ l(l992l) , seems to favour galaxies 
as the sites of production of the metals in the ICM. 

Accepting that the metals in the ICM originated in the 
E and SO galaxies of the cluster, two mechanisms exist to 
extract the newly synthesized elements from the individ- 
ual galaxies: the above mentioned GW a nd ram p r essure 
stripping. Arguments have been given bv iRenzinil l)l997h 
favouring the GW scenario, mostly based on the observa- 
tion that the "iron mass-to-luminosity ratio" is roughly 
constant independently of cluster richness and tempera- 
ture, while the ram pressure mechanism should be more 
efficient, extracting more metals for a given stellar con- 
tent, in richer clusters. 

Th ough the role of ram p ressure stripping is still de- 
bated ijMori fc Burkertll2000|) . from here on we will limit 
to the GW scenario for the pollution of the ICM, bearing 
in mind that the addition of other mechanisms (metal pro- 
duction in pre-galactic objects and/or ram pressure strip- 
ping) would allow to inject even more metals into the ICM, 
further favouring the enric hment. In our models, GWs are 
powered by SN feed-back l(Chiosill2000bh . as typically as- 
sumed in literature although AGNs have b een suggested 
as an alternative source of energy input IjRenzini et alJ 
Il993t iRomano et ahlEoO^ . 

Over the past two decades, modelling the metallicity 
of the ICM has been the subject of a great deal of stud- 
ies that we summarize briefly in the next section. In this 
paper we present a model for the chemical evolution of 
a cluster, in which galaxies eject part of their gas con- 
tent by the GW mechanism, thus playing the role of stars 
in the classical mod els for the chemical evolution of the 
interstellar medium IjTinslev 1980; Matteuccilll997h . 

The plan of the paper is as follows. Section 2 shortly 
reviews the problem of the gas and metal content of the 
ICM, and previous studies in literature. In Section 3 we 
suggest that a non-standard, variable initial mass function 
(IMF) for the stellar content of a galaxy could improve 



upon our understanding of the problem, and discuss the 
corresponding galactic models. Section 4 gives some simple 
estimates of the global properties of clusters, as expected 
from galactic models with the variable IMF or with the 
standard Salpeter IMF; we introduce the concept of intra- 
cluster mass-to-light ratio (ICMLR) as a measure of the 
amount of ICM gas. In Section 5 we introduce our new 
model for the chemical evolution of a cluster as a whole, 
the underlying analogy between the inter-stellar and intra- 
cluster medium, the concepts of galactic formation rate 
(GFR) and of galactic initial mass function (GIMF), and 
the model equations. Section 6 contains the detailed dis- 
cussion of a fiducial model, calibrated on the observed 
ICMLR and the present-day luminosity function (LF) of 
galaxies as the key constraint for the GFR and GIMF; 
a comparison between results with the Salpeter IMF and 
with the non-standard IMF is also made. In Section 7, a 
set of models is presented with different galaxy formation 
histories. In Section 8 we present cluster models computed 
on the base of galactic models especially selected to repro- 
duce the correct metal content of clusters. In Section 9 we 
discuss the predicted [a/Fe] abundance ratios in the ICM. 
Summary and conclusions are drawn in Section 10. 

2. Gas and metal in the ICM: previous studies 

Various early studies investigated whether "stan- 
dard" chemical models for galaxies ca n explain the 
amount of metals detected in the ICM Jvigrouxlll977t 
Iffinnes fc Biermannll98'ollDe YounJl978l) : by "standard" 
we mean a chemical model with the same physical ingredi- 
ents (mainly, stellar IMF and yields) suited to reproduce 
the Solar Neighborhood. Amidst these early studies, we 
menti on in particular the one by iMatteucci fc Vettolani 
(1988J) as the first attempt to link directly the metallic- 
ity of the ICM with the properties of the corresponding 
galaxy population. To this aim, the authors developed a 
modelling technique that has been widely adopted after- 
wards. Basing on a grid of models of elliptical galaxies 
with GW, they assigned to any given galaxy of final stel- 
lar mass M t , or equivalently of present-day luminosity 
L„, the corresponding masses of gas and iron ejected in 
the GW, Mgl s and Mj? e . Integrating these quantities over 
the observed LF, they calculated the total masses of gas 
and iron globally ejected by the galactic population in the 
cluster. Their main conclusions were: 

(1) the iron content of the ICM can be reproduced with a 
standard Salpeter IMF in the individual galaxies; 

(2) the global amount of gas ejected as GW is much 
smaller than the observed mass of the ICM, hence the 
ICM must be mostly primordial gas which was never in- 
volved in galaxy formation. 

Later on, this early success in reproducing M F c e M 
turned out to be favoured by the low gravitational poten- 
tial wells of model galaxies, calculated only on the base 
of their luminous, baryonic component. Once the poten- 
tial well of the dark matter halo is taken into account, the 
ejecta of SN la hardly escape the galaxy and the metal pol- 
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lution of the ICM by GWs is muc h reduced l|David et alJ 
1199114 iMatteucci fc Gibsonl 1 1995ft . In this case, standard 
chemical models fail to reproduce the metal content of the 
ICM. Some non-standard scenarios were thus invoked to 
solve the riddle, such as: 

- a more top-heavy IMF than the Salpeter one, 
with logarithmic slope x ~ 1.0 rather than the 
standard value x — 1.35 



Matteucci fc Gibson 




Il995l ICibson & Matteucci 


1997b a; 


lLoewenstein & Mushotzkvlll996ft: 



- a bimodal IMF with an early generation of massive 
stars heavily polluting the ICM, followed by a more 
normal st ar formation phase producing the sta rs we 
see today l)Arna,ud et a.lJh 992HElbay, et al.Hl 995h . 

These models, where SN II from massive stars play 
the main role in the metal pollution of the ICM, were 
further supported by the enhanced abundances of a- 
elements with respect to iron detected with ASCA 
l|Mushotzkv et al.lll996fr . 

Just as in the early work by IMatteucci fe Vettolani 
most authors conclude that GWs cannot account 
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for the huge amount of gas pre sent in the ICM, wh ich is 
2-5 times the mass in galaxies l)Arnaud et alJll992T) . The 
IC M must then con sist, 50 to 90%, of primordial gas. 

iTrenthaml l|l994lL on the base of the steep slope of the 
LF at the low luminosity end observed in clusters, sug- 
gested that all the intra-cluster gas could have originated 
in dwarf galaxies, since these are numerous in clusters 
and they are expected to eject a large fraction of their 
initial mass as GW, due to their s hallow potentia l wells . 
This s uggestion was discarded bv iNath fc Chibal (|l995h 
and bv iGibson fc Matteucci ( 1997bft . who calculated de- 
tailed models of dwarf galaxies and related GW ejection 
to show that galaxies cannot be the only source for the 
whole intra-cluster gas, even in the case of the steepest 
observed LF (hence the largest contributi on from dwarfs). 

Th e models were further refined by iMartinelli et al.l 
(2000J) who made use of one-zone and multi-zone GW 



models of elliptical galaxies and studied the dependence of 
the ICM abundances with redshift. While the abundance 
ratios [O/Fe] are in both cases within the observational 
uncertainties, the abundances [Fe/H] are very large and 
require large diluti ons by primo r dial g as. 

Very recently, IPipino et alJ ([2002) showed that the 
Salpeter IMF might reproduce the observed metallicity 
of the ICM, provided 100% efficiency of energy feed- 
back is adopted for SN la. In this case, all the iron pro- 
duced at late times escapes into the ICM in a continuous 
wind/outflow. The problem with this scenario, however, is 
that it inevitably predicts strongly sub-solar [a/Fe] ratios 
in the ICM. 

Finally, a recent attempt aimed at reproducing si- 
multaneously the iron abundance, th e ratio [O /Fe] , and 
the gas mass, was by IChiosi ll2000b|) who made use of 
multi-zone models of elliptical galaxies and adopted a non- 
standard IMF for their stellar content. More details on the 



models bv lChiosl (2000b) will be discussed in the next sec- 
tion as our present model stems from that work. 

3. A non-standard IMF 

As reviewed in the previous section, a wealth of work in 
literature suggests that some non-standard scenario (or 
IMF) must be invoked to account for the metals in the 
ICM. We recall that a non-standard IMF has been sug- 
gested for elliptical galaxies also on the base of other, in- 
dependent arguments: 

— a top-heavy IMF (x ~ 1.0) is better suited to 
reproduce the photome tric properties of ellipticals 
l|Arimoto fc Yoshii[l987ft : 

— systematic variations of the IMF in ellipticals of in- 
creasing mass might explain the increase of the mass- 
to-light (M/L) ratio with galactic luminosity, that 
is the s o-call e d "tilt of the Fundam e ntal Plane" 
dLarsor] Il98fit iRenzini fc Ciotti Il993t IZenf fc SillJ 
ll99fiHChiosi et alJll998l) . 

What physical reason may lead to a different IMF in dif- 
ferent situations? From the theoretical point of view, a 
turn-over of the IMF at low masses is expected, related 
to the Jeans scale (thermal support) and to the scale of 
magnetic support against gravit ational collapse (|LarsorJ 
Il998t IPadoan fc Nordlundl l20ol . The conditions of the 
ambient gas may thus influence the cut-off of the IMF at 
low masses, and hence the mass fraction of a stellar gen- 
eration that is locked into ever-living, very low mass and 
low luminosity stars. A variation of this "locked-up frac- 
tion" is crucial for the efficiency of the metal e nrichment 
pr oduced by a st ellar population l|Tinslevlll980l) . 

lLarsorJ ( 1998h suggested the following functional form 
of the IMF: 



dN 



d log M 



oc M x exp — — 



M 
~M 



where M s is a typical mass scale related to the Jeans mass. 
In brief, this IMF is a Salpeter power law down to a typical 
peak mass 

M p '- 

x 

below which there is an exponential cut-off. The peak mass 
varies with the temperature and density of the parent gas 
as expected from Jeans' law: 

M p ocT 2 

In warmer and/or less dense gas, therefore, the typical 
peak mass M r inc reases and the locked-up fra c tion is 
lower.lLarsonlljl998h . as well as lChiosi etHI l|l998l) : fchiosi 
(2000b), pointed out that a minimum temperature for the 
star forming clouds is set by the cosmic microwave back- 
ground, whose temperature increases with redshift, be- 
coming higher than the typical temperature of present- 
day molecular clouds at z > 2. At increasing redshift, the 
sole background then produces an increasing Jeans mass 
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and hence an IMF skewed toward massive stars. Other 
heating sources, like feed-back from massive stars or the 
UV background, and a reduced cooling rate at low metal- 
licities, may enhance the effect further. 

From the observational point of view, the issue of the 
variation of the IMF with ambient conditions is still open 
and very much debated. In local studies, some authors 
underline that, within the uncertainties, data are consis- 
tent with a constant IMF (Kroupa 2002), others find evi- 
dence of variation, e.g. between clust er and field llMasse 



19981) in clouds of d ifferent density l|Briceno et al.ll200 



Luhman et alJ^O Q j j), or betwee n the disk and the halo of 



the Milky Way <|Chabrierl l2f 

The IMF of the early star form ation activity at high 
redshi ft is even harder to probe. Irlernandez & Ferraral 
( 200l]) find indirect evidence in the local halo for a typical 
stellar mass scale (Larson-like) increasing with redshift. 
A similar scenario has been very recently advocated by 
iFinoguenov et ail l)2003h on the base of the different metal 
content in groups and clusters — and the different typical 
redshift of formation of their stellar content. Indications of 
a top-heavy IMF at high redshift (z > 3 — 6) has been als o 
found for Lyman Break Galaxies l|Ferguson et al.ir2002^) . 
All of this agrees with a number of recent theoretical 
studies suggesting that the first generat ions of stars were 
strongly skewed toward massive stars l|Abel et al.l l2002: 
Brommetal 1 l2002t IChiosll2000al and references therein). 

In view of these results and theoretical arguments it 
is certainly a legitimate working hypothesis to consider 
an IMF with a physical dependence on the environment, 
through the typical Jeans mass. 

3.1. Galactic models with the PNJ IMF 

A behaviour resembling the one de scribed above i s pre- 
dicted by the theoretical IMF by IPadoan et al ] ill997t 
hereinafter PNJ), which features a peak mass 



M p = 0.2 M Q 



(iok) 1 1000 < 



2.5 km s" 



where T, n, and a are the gas temperature, number den- 
sity, and velocity dispersion, respectively. Although the 
physical der ivation of the PN J IMF has been sometimes 
questioned l)Scalo et al.lll998|) . one can still adopt it in 
galactic models as a tentative recipe yiel ding the typica l 
behaviour expected for the Jeans mass; see lChiosll |2000b) 
for further discu ssion . 
IChiosi et all 



( 1998) developed chemo- 
thermodynamical models following the thermodynamical 
evolution of the gas in an elliptical galaxy, and the 
corresponding variations in the IMF according to the 
PNJ recipe. The characteristics and behaviour of these 
models as a function of galactic mass and redshift of 
forma tion a re disc u ssed in full details in IChiosi et alJ 
l|l998|) and IChiosil l)2000b|) . Here, we briefly underline 
the qualitative trends with respect to (a) mass and (b) 
redshift of formation. 



(a) At increasing galactic mass, the average density of the 
object p decreases and the typical peak mass M p in- 
creases, yielding a lower locked-up fraction. 

(b) At increasing redshift of formation Zf or , the tem- 
perature of the protogalactic gas increases, since it 
can never fall below the corresponding temperature of 
the cosmic microwave background, T > TcMB{zfor), 
which increases with redshift; hence, the peak mass M p 
is higher and the locked-up fraction is lower. 

We remark here that the above mentioned trends are by 
no means drastic: the peak mass M p hardly exceeds 1 M , 
and the "high M p " phase is limited to the early galactic 
ages; after the initial stage, in fact, the system reaches 
a sort of thermodynamical balance, with the peak mass 
and the IMF settling on quite standard values. The over- 
all pictur e loosely resembles the bimodal behaviour sug- 
gested by Elbaz et al . (1995). with an early phase domi- 
nated by massive stars followed by a more normal star for- 
mation phase producing the low-mass stars we see today. 
However, in our models the IMF naturally and smoothly 
varies in time following a physical prescription, rather 
than an imposed bimodal behaviour. 

Though not long-lasting, the variations in the early 
phases suffice to differentiate the resulting galactic models, 
making them successful at r eproducing many features of 
observed ellipticals, such as l|Chiosi et al.lll998|) : 

— the tilt of the Fundamental Plane; 

— the analogous of the "G-dwarf" problem, or the lack 
of a large population of low metallicity stars, de- 
tected in the spectral energy distribution of ellipticals 
(|Bressan et al.lll994l IWorthev et al.lll996l): 

— the high fraction of white dwarfs ijBica et alJ ll996): 

— both the colour-magnitude relation and the trends in 
a-enhancement with mass at the same time, thereby 
overcoming the well-known dichoto my between the 
"classic" an d "inverse" GW scenario l Matteuccil fl992. 

Il994lll997l) . 

This last point is worth commenting further, as the mod- 
elling of the GW influences directly the predictions con- 
cerning the ICM. GW models of elliptical galaxies with a 
constant IMF (whether Salpeter or more top-heavy) face 
the following puzzle. The colour-magnitude relation sug- 
gests that GWs occur later in more massive ellipticals than 
in smaller ones, so that star formation and chemical en- 
richment proceed longer and the stellar population reaches 
redder colours in more luminous objects. On the other 
hand, metallicity indices, if interpreted as abundance in- 
dicators, suggest that the [Mg/Fe] ratio increases with 
galactic mass; this requires GWs to occur earlier in more 
massive galaxies, where SN II should dominate the chem- 
ical enrichment to make the resulting abundance ratios in 
stars a-enhanced. 

This dichotomy between the so-called "classic" and 
"inverse" GW scenario, hampers predictions of the metal 
pollution of the ICM, since two competing sets of GW 
models are to be considered. It is therefore quite appeal- 
ing, when we address the chemical enrichment of the ICM, 
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that the variable IMF scenario described above can repro- 
duce both observational constraints, with a unique set of 
models. 

IChiosil ll2000hl) first analyzed the predicted metal pol- 
lution of the ICM when galaxy models with the PNJ IMF 
are adopted. To this purpose, he calculated multi-zone 
chemical models of ellipticals with the PNJ IMF. 

The adoption of radial multi-zone models, rather than 
simple one-zone models, has in fact important conse- 
quences on the predicted enrichment o f the ICM, as un- 
derlined first bv lMartinelli et all <HpOO). When the radial 
structure of an elliptical galaxy is considered, with the cor- 
responding gradients in density, colours etc., it turns out 
that the GW does not develop instantly over the whole 
galaxy, but it tends to set in earlier in the outskirts, where 
the potential well is shallower, and later in the central 
parts. This means that star formation and chemical en- 
richment proceed longer in t he centre than in the outer 
regions l|Tantalo et a.l.lll99St iMartinelli et al~lll998^ . and 
the GW ejected from different galactic regions is metal 
enriched to different degrees. 

The models calculated bv lChiosil l)2000tj) account for 
this effect by dividing the galaxy into three zones: a central 
sphere where star formation and metal production is most 
efficient and lasts longer; an intermediate shell where the 
GW sets in earlier and the metal production proceeds to 
a lesser extent; an outer corona where the gas is expelled 
almost immediately, with virtually no star formation and 
chemical processing. This behaviour is the combined result 
of the shallower potential well when moving outward in 
the galaxy (as in standard models with a constant IMF) 
and of the varying M p in t he PNJ IMF w hen moving to 
outer, less dense regions; see lChiosl (|2000bl) for a detailed 
discussion. 

For the sake of comparison, analogous models with 
the Salpeter IMF were also calculated, with mass range 
[0.18—1201 Mn, as suited to model elliptical galaxies 
l|Tantalo et alJll99S|) . 

In all these models l|Chiosi et al.lll998t IChiosill2000r^ 

the metal production and recycling, and the resulting 
abundances of the GW ejecta, are follo wed with the 
chemi cal evolution network developed bv IPortinari et all 
l)l998|) . Chemical yields of massive stars were derived 
from the stellar tracks of the Padua group for the pre- 
supern ova phases, and th e n link ed to the SN II mod- 
els by IWooslev fc Weaver! l|l995l) . rescaled to the same 
core masses. Yi elds for low and intermed iate mass stars 
were taken from lMarigo et all l|l996tll998(l . however these 
are of minor importance in the present work: for the 
GW ejecta and the enrichment of the ICM, most impor- 
tant are iron and a-elements, produced by supernovae. 
The chemical network also includes type la SN, which 
are important iron contr ibutors, adopting the recipe by 
iGreggio k. Renzinl Jl983l) for the rate, and t he ejecta from 
the W7 model bv iThielemann et ail jl993|) . We refer the 
reader to the original papers for further details on the 
chemical network. 
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Fig. 1. Mass fractions in GW and in remaining galaxy 
as a function of the initial (proto)galactic mass M, for 
four different redshifts of formation. Thick lines: galactic 
models with the PNJ IMF; thin lines: galactic models with 
the Salpeter IMF. 
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Fig. 2. Iron abundance (in solar units) of the GW, for 
galaxies of initial (baryonic) mass M as indicated in each 
panel, and as a function of the redshift of formation. Thick 
lines: galactic models with the PNJ IMF; thin lines: mod- 
els with the Salpeter IMF (redshift independent). 
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3.2. Galactic ejecta: PNJ vs. Salpeter IMF 

For a better understanding of the results concerning the 
chemical evolution of the ICM, let's first inspect the pre- 
dicted GW ejecta of the galactic models when the variable 
IMF or the Salpeter IMF are adopted in turn. 

A (proto)galaxy of initial baryonic mass M formed at 
redshift Zf or ejects a mass Eqw{M, Zf or ) of gas as GW, 
while a mass 

R ga l(M, Zf or ) = M - E GW (M, Zf or ) 

remains as the baryonic component of the galaxy we "see 
today", that is as the stars (and remnants) of the final 
galaxy. For galactic models with the Salpeter IMF, the 
various quantities depend only on M and not on Zf ori as 
there are no temperature effects on the IMF in that case. 
For the purpose of the chemical evolution of the ICM, 
Rgai can be viewed as a "galactic remnant" , passively sub- 
tracting mass from further chemical/galactic processing, 
analogous to stellar remnants in chemical models of indi- 
vidual galaxies. The gas shed by long-lived stars after the 
GW episode is supposed to remain in the galaxy and not 
contribute further to the enrichment of the ICM. 

In Fig. ^ we plot the mass fraction of gas ejected in 
the GW 

, Eqw 

and the complementary mass fraction locked into the 
"galactic remnant" , 

t Rgal -. t 

J gal — = 

for galactic models with the variable PNJ IMF and for 
models with the Salpeter IMF — thick and thin lines, 
respectively. Mass fractions refer to the total initial bary- 
onic mass of the (proto)galaxy. The amount of ejected gas 
is larger in the case of the PNJ IMF, since less mass is 
locked into low-mass stars thanks to the high M p in the 
early galactic phases. In particular, with the PNJ IMF the 
mass ejected as GW is always larger than that locked up 
in the galaxy. The difference with the Salpeter case gets 
sharper for larger (proto)galactic masses, and for higher 
redshifts of formation (as already mentioned, models with 
the Salpeter IMF bear in fact no dependence on the red- 
shift of formation). It is worth underlining here the follow- 
ing "inverse" behaviour of the models with the PNJ IMF 
with respect to what is generally found from models with a 
constant IMF. According to the general consensus, larger 
galaxies store a larger fraction of their mass into stars 
and eject a lower fraction of gas in the GW, while smaller 
galaxies, due to their shallower potential wells, are more 
efficient in wind ejection. This trend is indeed evident in 
the Salpeter galactic models of Fig. ^ The models with 
the PNJ IMF, on the other hand, show quite the opposite 
behaviour: larger galaxies eject a larger fraction of their 
initial mass in the wind and lock a lower fraction into stars, 
due to the higher peak mass that characterizes them in the 
early phases. This IMF effect overwhelms their deeper po- 
tential wells, and the trend becomes stronger and stronger 



with increasing redshift of formation. There are in fact ar- 
guments advocating large baryon losses from galaxies in 
general (ISilkll2002fl . 

The masses of iron and oxygen, respectively, ejected 
in the GW are Egw,Fc and Eqw,Oi so that the metal 
abundances of the GW are given by: 

v Ecw,Fe v E G W,0 

-&GW.Fe — — ^ -X-GW,0 = — ^ 

J^GW J^GW 

The quantities Eqw, Egw.Fc, Egw.o for a grid of values 
of (M, Zfor) are tabulated in Chiosi (2000, his Table 4, 
detailing the contribution of the three different shells). 

Fig. |3 shows the iron abundance X G w,Fe m the gas 
ejected as GW, again comparing the Salpeter IMF (thin 
lines) and the PNJ IMF case (thick lines). In most cases, 
the galactic ejecta in the PNJ models are more metal-rich 
than in the Salpeter case, up to a factor of five or more 
in the case of the most massive galaxies, and for high 
redshifts of formation. In the PNJ models, in fact, more 
galactic gas gets recycled through massive stars, effective 
metal contributors, and less mass gets locked into low- 
mass stars, before the GW occurs. 

From the trends described above, we expect that mod- 
els of ellipticals with the PNJ IMF predict a more efficient 
metal enrichment of the ICM, and a higher fraction of 
its gas originating from GWs, with respect to "standard" 
models . The fi rst results in this respect were discussed by 
IChiosil ll2000bh . 

3.3. M/L ratio and IMLR of galaxies: PNJ vs. Salpeter 

A popular way of measuring the efficiency of metal pro- 
duction of cluster galaxies is the so -called iron mass t o 
light ratio (IMLR) introduced by ICiotti et all l)l99l|) : 
iRenzini et alJ l|l993(L This quantity is ana logous to the 
"global yield " in chemical evolution models l(Tinslevl ll980: 
|Pa,gellll997t) . defined as the ratio of the global amount 
of metals produced by a generation of stars to the mass 
that remains locked in remnants plus low-mass, ever liv- 
ing stars. In the case of clusters, the mass locked in stars, 
i.e. in the stellar component of galaxies, is replaced by 
their luminosity (in the B-band), which is directly mea- 
surable, while the estimate of the mass would be indirect 
and would require some as s umpt i on ab out the M/L ratio 
of cluster galaxies. iRenzinil l)l997l l200 0i) estimates, for the 
IMLR of the ICM: 

M Fe i 
= (0.02 ± 0.01) h~* 

Lb 

where M[q M is the mass of iron in the ICM gas, Lb is the 
global luminosity of cluster galaxies, and ft.50 is the Hubble 
constant in units of 50 km sec -1 Mpc -1 . The I MLR ap- 
pears to be quite constant among rich clusters l)Renzinil 
119971) . More recent estimates of the IMLR taking into ac- 
count the existence of radial gradients of iron abundance 
in the ICM and considering gas masses, iron masses and 
luminosities within the same radius, have slightly lowered 
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the estimate l|Finoguenov et alJl2000tl200l|) : 



1V1 ICM 



(0.01 -0.015) ft 



50 



(1) 



These latter estimates have the advantage to be car- 
ried out consistently over the same cluster volume 
Metallicities are in fact typically m easured at most 
within half of the virial radius (Fino guenov et al 



2000l l200lllDe Grandi fc Molendl200lUlrwin fc Bregmar 
2001); measurements hardly reach t^q, which is typically 
63% of r?nn . and in terms of cluster mass M 500 ~ 0.63A/ 2 oo 
as well (see iReiprich fc Bohringerl l2002'). 1 For the sake of 
the chemical enrichment what matters is the mass aver- 
aged metallicity, obtained by convolving metallicity dis- 
tributions with gas profiles within the same radius, and 
the IMLR should also be evaluated from the gas, met- 
als and galactic lumino s ity within t he same radius (as in 
iFinoguenov et alJ l2000l l200ll 12003^ . Since the observed 
metallicities probe only half of the cluster mass, consid- 
ering them representative of the whole gas mass (derived 
from the gas profile extended out to the virial or Abcll 
radius) is a significant extrapolation, especially in view of 
recent results on gradients: relaxed clusters with cooling 
flows systematically show iron abundance gradients (also 
besides the central peak) and even the metallicity distribu- 
tion in non-cooling flow clusters, though consistent with 
being uniform, is better fitt ed with a negative gradient 
llDe Grandi fc MolendilEool . 

The IMLR could also be overestimated in the presence 
of a substantial intra-cluster stellar populations. Most 
studies limit this diffuse population to a 10-20% of the 
total stellar cont ent, but values up to 4 0% have been re- 
cently suggested l|Arnaboldi et al JI2Q0J. 

Finally, we remind that red or NIR bands are a better 
probe of the actual star mass in galaxies, for they are less 
sensitive to recent sporadic star formation and more to 
the old underlying population. Hence the IMLR would 
be a better probe of the real "yield" of galaxies, if it were 
expressed in terms of the R to K band luminosity. This is 
becoming possible nowadays, as cluster LFs in red or IR 
band s are presently becoming available llMobasher et all 
2003t iDriver et all Il99i iMobasher fc TrenthanT 1998 : 
Trentham fc Mobasherl 1 19981 foe Propris et alJ 199S : 
Andreon fc Pelloll2000t lAndreon fc Cuillandrdl2002|) : see 
m fact iLin et all faOQ 1 



While these developments are certainly interesting for 
the future, here we will consider as our observational con- 
straint the "canonical" value of the IMLR in the B band 
given by QJ. 

Before modelling the cluster and its IMLR as a whole, 
it is worth considering what is the IMLR of the galactic 



1 We indicate, as customary, with rjv the radius correspond- 
ing to an overdensity iV times larger than the average back- 
ground density; r2oo is usually identified with the virial radius. 
Mm indicates the mass enclosed within rjv. 
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Fig. 3. The M/L ratio in the B band for living stars in the 
PNJ galaxies {dots and triangles, two example models, see 
legend) and in the Salpeter galaxies {thin line; mass limits 
for the Salpeter IMF are [0.18-120] M Q ). 



wind of individual model galaxies: 



IMLR GW 



E, 



GW,Fe 

L B 



(2) 



comparing the PNJ IMF to the Salpeter IMF case. The 
GW ejecta Ecw,Fe have been discussed above, now we 
discuss the B-band luminosities Lb of our galaxies. A 
galaxy, or "galactic remnant" R ga i left over after the galac- 
tic wind, consists of three baryonic components: the living 
stars producing the galactic luminosity, the dark stellar 
remnants and the gas shed by stars after the GW. The 
first component progressively loses mass in favour of the 
other two; namely, within R ga i the fraction F* of living 
stars decreases with time (Fig. 0}. 

Comparing the PNJ and the Salpeter models, two con- 
trasting effects contribute to determine the final luminos- 
ity of a galaxy. With the PNJ IMF, in the early galactic 
phases the typical stellar mass was skewed to higher val- 
ues and less material was locked into low mass, very low 
luminosity stars. As a consequence, the living stars are 
on average more massive and more luminous in the PNJ 
models than in the Salpeter models, so that their typi- 
cal M/L ratio is lower (Fig. |3J). On the other hand, with 
the PNJ IMF more mass went into remnants in the early, 
top-heavy star formation phases, and the gas restitution 
fraction is higher (even after the GW) because of the lower 
number of ever-living low-mass stars; as a consequence, 
living stars often represent a lower mass fraction of the fi- 
nal galaxy in the PNJ models, especially for high redshifts 
of formation and/or large masses (Fig.QJ. Globally, the to- 
tal M/L ratio of the galaxies, or "galactic remnants" made 
of stars+remnants+gas (ejected by stars after the GW), 
is shown in Fig. |3J With respect to the Salpeter models, 
models with the PNJ IMF show a composite behaviour: 
at low redshifts of formation {zj or < 10) they have lower 
M/L ratios, for the effect of having more luminous stars 
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Fig. 4. Fraction of living stars within the "galactic rem- 
nant" Rgai , decreasing with age in favour of remnants and 
ejected gas. PNJ models for three initial (proto)galactic 
masses and two redshifts of formation shown for the sake 
of example (see legend). Thick lines with no dots corre- 
spond to the Salpeter models. 




age [Gyr] 

Fig. 5. M/L ratio of the global "galactic remnant" R ga i 
(stars+remnants+gas ejected after the GW) , as a function 
of age. Symbols as in Fig. 0] 



prevails; at high redshifts of formation, they have higher 
M/L ratios, for the effect of the larger amount of remnants 
prevails. 

Finally, Fig. displays the IMLRgw of the galax- 
ies (Eq. [2J as a function of the initial mass M of the 
(proto)galaxy and of its redshift of formation Zf or ] thick 
lines are for galaxies with the PNJ IMF, thin lines are for 
the Salpeter case. Notice that IMLRgw depends on Zf or 
also for the Salpeter case, for although Eaw,Fe is indepen- 
dent of redshift with the Salpeter IMF, Lb decreases due 
to increasing galactic age at increasing Zf or (the effect is 
however negligible between Zf or — 20 and 5). 
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Fig. 6. The IMLR in the GW for model galaxies with the 
PNJ IMF (thick lines) and the Salpeter IMF (thin lines). 



The trend of IMLRgw with mass M is opposite in 
the two cases, reflecting the effects already commented 
upon with Fig. 1 and 2. At small masses, models with the 
Salpeter IMF eject more gas and metals in the GW, retain- 
ing less mass in stars and having a lower final luminosity; 
their IMLR GV k is thus higher than that of the PNJ IMF 
models. At large masses, say M > a few 10 10 Mq, the be- 
haviour is reversed and it is the galactic models with the 
PNJ IMF that eject a higher amount of metals relative 
to the luminous mass retained in the "galactic remnant" . 
For Zf or < 5, GW ejecta are the same as in the Zf or = 5 
the background temperature has by then dropped 
to low values and there is no further evolution on its basic 
effects on the IMF. Therefore for z for < 5, for both IMFs 
the corresponding IMLRgw changes (decreases) just be- 
cause of the effects of increasing Lb at younger ages. 

Most important for the sake of the chemical enrich- 
ment of the ICM, notice that the Salpeter models for 
all masses have an intrinsic IMLRgvk < 0.01, i.e. always 
lower than what is observed in clusters (Eq. flj . We can 
thus anticipated that model galaxies with the Salpeter 
IMF are incapable of reproducing the high IMLR observed 
in clusters. On the contrary, models with the PNJ IMF 
at high masses and/or for high redshifts of formation do 
reach IMLR GW > 0.01. 

To predict the IMLR in the cluster as a whole, we need 
to convolve the respective IMLRciy with the Galactic 
Formation History of the cluster; we will see (|0J| that 
models with the PNJ IMF yield a higher global IMLR, 
in better agreement with observations. This is due to the 
fact that small galaxies, although dominating the LF in 
number, actually represent a minor contribution in mass, 
so that the main role in the ICM po llution is playe d by 
galaxies with M > 10 10 M Q (see also lThomaslll999|) . For 
the latter type, galactic models with the PNJ IMF have 
in fact an IMLRgh/ of the order of 0.01-0.02, comparable 
to the typical value ^0.01 for clusters. 
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4. Global cluster results: the monolithic approach 

As mentioned in ^2 the most popular way to calcu- 
late the metal content of the ICM of a cluster on the 
base of its population of elliptical galaxies, is to de- 
velop a grid of galactic models and integrate their GW 
ejecta over the observed LF, a method first introduced by 
iMatteucci fc Vettolanil 1 19881) . In this approach, all the el- 
lipticals in the cluster are assumed to be coeval, hence we 
will call this the "monolithic" approach. 

In this section we adopt the monolithic method, with 
the PNJ and the Salpeter galactic mode ls in turn. Our ref - 
erence LF is the observed B-band LF bv lTrenthaml l|l998l) , 
since it is very deep (down to magnitude Mb = — 11) and 
since the observed IMLR is also referred to the B-band 
luminosity (see <J3 .3|> . This LF, displayed in Fig. is a 
weighted mean of the LFs of 9 clusters at low redshift 
[z < 0.2). Besides showing the standard exponential cut- 
o ff at brigh t magn itudes (Mb < —20), characteristic of 
a ISchechterl l|l97fil) function, this LF also steepens at the 
faint end (Mb > —14), that is in the regime of dwarf 
galaxies. 

The luminosity evolution of the "galactic remnants" is 
known as a function of the initial (proto)galactic mass M 
and of the redshift of formation Zf or (SJ22J- As to the 
relation between the redshift Zf or and the correspond- 
ing age of the galaxy, we adopt a flat A Cold Dark 
Matter (ACDM) cosmology with 0^=0.3, fi A =0.7 and 
Hq = 65 km sec -1 Mpc -1 . Once a value for Zf or is fixed, 
each luminosity bin of the LF corresponds to some initial 
galactic mass M, and we sum the contributions of the dif- 
ferent bins to estimate the global amounts of ejected gas 
and metals. 

The GW ejecta of the PNJ models are quite sensitive 
to the exact epoch (redshift) of formation of the individual 
galaxies, as discussed in E I3.2I and we perform the exercise 
for Zf or — 5, 10, 13, 15, 20. (The case z/ or =13 is obtained 
by interpolation in the grid of galactic models). For the 
Salpeter IMF, instead, computing the case Zf or =5 is suf- 
ficient because the GW ejecta are fixed, and the age and 
luminosity differences for Zf or >5 are negligible (Fig. EJ. 

Results for global integrated quantities with the PNJ 
and the Salpeter models are listed in Tabled The 3 rd col- 
umn shows the global IMLR for the cluster. The Salpeter 
models fail to reproduce the observed IMLR by an order 
of magnitude or so. The PNJ models perform much bet- 
ter, with an iron production and hence an IMLR increas- 
ing with redshift of formation, as expected from H3.3I For 
z/ or =10-13 the IMLR falls in the observed range (Eq.[TJ, 
while models with Zf or >13 have too large an iron pro- 
duction with respect to observations. 

4.1. The contribution of Type la supernovas 

At increasing redshift of formation, the PNJ IMF is more 
and more skewed toward massive stars in the early phases, 
implying a lower and lower ratio of the number of SN la 
vs. SN II. Column 8 in Table ^ reports the fraction of 
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Fig. 7. The observed B-band luminosity function of clus- 
ter galaxies by Trentham (1998), in relative frequency. 



the global iron produced and ejected into the ICM due to 
SN la; as expected, it decreases with increasing redshift, 
from 55 to 35%. The Salpeter IMF has a larger contribu- 
tion (65%) from SN la than any case of the PNJ IMF. 

The same trends are reflected in the [O/Fe] ratio of 
the ejected gas (9 column in the Table), increasing with 
redshift. Within the present uncertainties about the [a/Fe] 
ratio in the ICM (©, all of the values for [O /Fe] in Table[T] 
are acceptable. For the "favoured" cases (on the base of 
the IMLR) of z/ or =10-13, the resulting [O/Fe] ratio is 
close to solar. 

One concern when computing the iron production and 
ejection from galactic models is the uncertainty in the rate 
of SN la, due to the poorly known evolution of the pro- 
genitors. As a consequence, one may wonder if the high 
iron abundances in the ICM can be reproduced by adopt- 
ing a higher rate of SN la rather than a top-heavy or 
non-standard IMF. However, the [O/Fe] ratios we obtain 
indicate that the relative contribution of SN la vs. SN II 
predicted for the ICM are grossly correct. 

With the Salpeter IMF and standard GWs models of 
elliptical galaxies, the predicted contribution to the iron 
enrichment of the ICM is short by 5-10 times with respect 
to the observed IMLR. If we were to compensate for this 
discrepancy by adjusting the SN la rates to the required 
level — or by assuming that later iron production from 
SN la escapes the galaxy after the main GW episode - 
the corresponding [O/Fe] ratios would decrease by 0.7-1 
dex, by far too low with respect to observations. 

4.2. The problem of dilution 

Besides the amount of ejected metals as traced by the 
IMLR, an important question concerns the global gas mass 
ejected from galaxies as GW: can this account for the 
whole of the ICM, or conversely most of the ICM is primor- 
dial gas that was never involved in galaxy formation? The 
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IMF 


Z f or 


IMLR 


M gal 

Ltot 


M G w 
Ltot 


M GW 
M gal 


X Fc ,GW 
X Fe,® 


fsNIa 


[O/Fe] 


M IC M 
M G w 


Mj CM 
M gat 


X Fc ,ICM 


PNJ 


5 


0.007 


7.2 


11.1 


1.54 


0.358 


0.53 


-0.08 


3.3 


5.1 


0.11 


PNJ 


10 


0.009 


8.2 


13.0 


1.59 


0.376 


0.54 


-0.07 


2.8 


4.5 


0.13 


PNJ 


13 


0.018 


10.6 


21.1 


1.99 


0.477 


0.47 


+0.07 


1.8 


3.5 


0.27 


PNJ 


15 


0.029 


12.5 


28.9 


2.31 


0.562 


0.43 


+0.16 


1.3 


2.9 


0.44 


PNJ 


20 


0.065 


17.4 


50.6 


2.91 


0.714 


0.36 


+0.21 








Salpeter 


5 


0.002 


8.6 


6.2 


0.72 


0.143 


0.65 


-0.09 


6.0 


4.3 


0.02 



Table 1. Results for integrated quantities for the cluster, in the monolithic approach. 1 st column: IMF of model 
galaxies. 2 nd column: redshift of formation of galaxies. 3 rd column: global IMLR for the cluster. 4 column: global 
M/L ratio for cluster galaxies (ratio between the global mass in "galactic remnants" and the global B-luminosity) . 
5 th column: ratio between the mass of gas ejected as GW and the global B-luminosity of cluster galaxies. 6 th column: 
ratio between the mass in gas ejected as GW and the mass in "galactic remnants" . 7 th column: iron abundance with 
respect to solar in the gas ejected as GW. 8 th column: fraction of the global iron production contributed by SN la. 
9 th column: [O/Fe] ratio in the gas ejected as GW (and hence in the ICM). 10 th column: "dilution factor" necessary 
to recover the observed ICMLR=M/cM/-ktot=37 MQ/L & (see text). 11 th column: ratio between the total mass of ICM 
gas and the mass in galaxies. 12 th column: average iron abundance in the ICM gas. 



observational constraint here is the global amount of ICM 
gas vs. the mass in galaxies; unfortunately, this quantity 
is relatively poorly determined. In literatures, values can 
be found in the range Micm /M ga i=2-20; a closer inspec- 
tion reveals that this large range can be mostly imputed 
to different assumptions for the M/L ratios for galaxies. 
In fact, while the mass in the ICM can be estimated di- 
rectly from X-ray observations, the mass in galaxies is 
estimated from their global luminosity, usually in the B- 
band, and assuming a M/L ratio. We remark here that for 
the chemical enrichment of the ICM, what is meaningful 
is the baryonic mass in galaxies, not their global mass or 
M/L ratio inclusive of the dark halo component. The real 
observable is Micm / Lb, arid with respect to this quantity 
results in literature a re much more homogeneous. 

From the study of iDavid et al of Hydra A one 

_ i 

finds Mjcm/Lv = 33 /i 50 2 Mq/L q ; the typical colours of 
ellipticals correspond to Ly ~ 1.3 Lb and the above entry 
translates into: 



Micm 



Me 
Lq 



lArnaud' et all lll992t assume Ly ~ 1.2 Lb for the typical 

colours of E and SO galaxies, and consequently derive for 

_ i 

their clusters Micm / 'Lv : (e+so) ~ (20 — 50) h 50 2 , with an 

average of 36 h 50 2 (see also lElbaz et al.ll995h . correspond- 
ing again to: 

< Mwm_ > _ 3Q/ ^i Mq w . th arange ( 17 _ 42 ) 
Lb Lq 



The analysis of the Coma cluster bv lWhite et al ] lll993h 
yields as well: 

Lb L & 



ICirimele et all ( 1997|) estimate Micm /M ga i for their sam- 
ple under the assumption of a typical M/Ly = 8 Mq/Lq 



for elliptical galaxies; their tabulated values correspond 
to Micm / Ly in the range 8.5-37.6 with an average value 
of 21 Mq/Lq. Taking again a typical Ly ~ 1.3 Lb for 
ellipticals, and considering that Cirimele et al. assumed 
Hq = 50, their data imply: 

<^>~i9/r^ 

Lb Lq 



In the sample bv iRoussel et al.l l)200fl) . the average over 
all the objects is: 

<^>~3i/r^ 

Lb Lq 

while for the hot clusters only, a somewhat higher value 
is obtained: 

Lb Lq 

although it depends on the volume probed: 
Micm/L b =27, 35 or 44 h~* Mq/Lq, for r < r 2 ooa,r 5 oa 
or r < r2oo respectively, since the gas distribution is 
ge nerally more extended tha n the galaxy distribution. 

iFinoguenov et alJ l|2003l) indicate a typical value of 
Micm /Lb ^30 for hot clusters within 0.4 rioo, adopting 
#o=70, or: 

Micm 



L, 



In summary, data indicate a typical value of "intra-cluster 
mass to light ratio" 

„ Mjcm i Mq 

ICMLR= ° - 30/1^2 

Lb Lq 

with somewh a t lower values in th e sam ples by 
ICirimele et ail l)l997|) : IFinoguenov et alJ l|2003|) (though 
still growing with radius in the latter) and somewha t 
higher values in the sample by IRoussel et al] l|2000|) . 
Whether the fraction of hot vs. cold baryons, that is the 
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efficiency of galaxy formation, increases noticeably with 
cluster temperature and richness is on the other hand 
still a matter of debate llBrvanll200(i iBalogh et al.ll200lt 
ILin et aklEool) . 

We favour MjcmI^b as a constraint rather than 
Micm /Mg a i, since the former is directly measured. Just 
as the IMLR is supposedly a more straightforward esti- 
mate of the efficiency of metal enrichment than the "clas- 
sical" yield M p e /M* (see ffOty . the ICMLR is a more di- 
rect constraint than the "classical" gas fraction Micm /M* 
because, like the IMLR, it is independent of a priori as- 
sumptions on the M/L ratio of galaxies. For our present 
favoured value of ft. = 0.65, we regard ICMLR^37 M@ / Lq 
as our constraint for the total amount of ICM gas in clus- 
ters. 

In the last three columns of Tablcfjwe give the dilution 
factors (i.e. the ratio between the total ICM gas mass and 
the mass provided by the GWs) necessary to recover the 
observed ICMLR, the corresponding ratio between ICM 
mass and mass in galaxies (the latter inclusive of stars, 
remnants and gas shed by stars after the GW) and the 
average metallicity in the global ICM. Typical dilution 
factors are in the range 2-3 for the PNJ models, and the 
ICM mass is typically 3-5 times the mass in galaxies. 

In the case of Zf or — 13, that is the case with the best 
results for the IMLR, roughly half of the ICM is expected 
to come from GWs, its average metallicity is close to the 
observed characteristic value of 0.3 solar, and the ICM 
mass is 3.5 times the mass in galaxies. Galactic models 
with a redshift of formation Zf or > 15 are ruled out since 
they eject not only too much metals (as marked by the 
high IMLR) but also too much gas in the GW to be com- 
patible with the observed ICMLR. We notice also that 
Salpeter models require a high degree of dilution (a factor 
of 6) and the corresponding metallicity in the ICM gas is a 
factor of 10 too low, as expected from the correspondingly 
low IMLR. 

Our global M/L ratio for galaxies (including all the 
baryonic components, living stars, remnants and gas shed 
by stars at late times) are in the range 7-12 for acceptable 
models with Zf or < 15, and ~8 for the Salpeter models. 
Correspondingly, the observed ICMLR=37 M /L im- 
plies that the mass in the ICM is 3-5 times the mass in 
galaxies. This is a factor of 2- 3 lower than t he w idely 
quoted value of 10 derived by I White et~al ] lll993h for 
h = 0.65. The difference just stems from their much lower 
M/Ls=6.4 h, adopted on the basis of dynamical argu- 
ments; this is lower than what is expected for the typi- 
cal stellar population in an old elliptical galaxy (cf. the 
M/L ratio of Salpeter models). The difference is quite ir - 
relevant for the problem discussed bv I White etall |l993), 
that is the baryon fraction in clusters: the baryonic mass 
is dominated by the hot ICM gas mass anyways, and is 
not much affected by uncertainties in the mass in galax- 
ies — though the latter might not be as negligible, as 
sometimes assumed. However, for the sake of the chem- 
ical enrichment of the ICM the effect is quite crucial. 
Adopting the ICM-to-galaxy mass ratio bv IWhite et alJ 



lll993h . E enzm 3 <l!997l) estimates the ratio between the 
metals in the ICM and that locked in the stars of cluster 
galaxies to be 1.65 fr~ 3 / 2 , or - 3 for h = 0.65. If the M/L 
ratio in galaxies is higher, as required e.g. by stellar popu- 
lation models for ellipticals, the "metal balance" is much 
less skewed toward the ICM so that the amount of metals 
in the ICM becomes comparable to that in galaxies. This 
underlines the importance, for a consistent modelling of 
the chemical evolution of the ICM, to adopt the observed 
ICMLR as a constraint rather than some independently 
derived ICM-to-galaxy ratio relying on external assump- 
tions on galactic M/L ratios. 

Finally, we remark that the actual ICM-to-galaxy 
mass ratio is important for the sake of explaining the 
"entropy floor" in low-T clusters, whether it requires 
strong supernova pre-heating or whether it is partly due 
to the removal of low-entropy gas by galaxy formatio n 
l|Brvanll2000t IBalogh et afll200lUTornatore et alJl2003l) 
Intere stingly, the recent cluster simulations bv lVaMaTnml 
(2002J) with star formation and self-consistent chemical 
evolution, suggest that the metallicities of the ICM can be 
reproduced provided the IMF is top-heavy with respect 
to Salpeter, and the corresponding mass in cold baryons 
(galaxies, or "star particles" in the simulations) is large 
enoug h to favour th e scenario of removal of low-entropy 
gas hv lBrvaril <|2000h. 

5. A chemical model for the ICM evolution 

The monolithic approach of ^assumes all the ellipticals in 
the cluster to be coeval, with an age of 13-14 Gyrs. This is 
a somewhat naive approach in the light of current under- 
standing of structure formation in the Universe. In the hi- 
erarchical scenario the typical size of virialized objects in- 
creases at decreasing redshift. There are various reasons to 
argue that this scenario, developed for the gravitationally 
dominating, collisionless dark matter, might not hold for 
baryons, whose evolution at galactic scales is probably de- 
coupled from that of dark matter — the angular momen- 
tum problem of disc galaxies in cosmological simulations, 
the observational evidence of early formation of ellipti- 
cal galaxies, the colour-magnitude relation of both disc 
and elliptical galaxies suggesting the stellar populations 
in larger objects to be older at odds with expectations 
from the hierarchical bottom-up assembly, etc. Hence it 
is possible that even in the hierarchical scenario, an old 
and roughly coeval stellar population of ellipticals can 
be accommodated, by suitable account of detailed baryon 
physics. Nevertheless, we explore now a picture in which 
the formation of cluster galaxies of different masses is red- 
shift dependent; as the extreme alternative to the mono- 
lithic approach, we will assume that a galaxy forms at the 
epoch when the corresponding dark halo is predicted to 
virialize according to the Press-Schechter formalism. We 
develop a chemical model following the epoch of forma- 
tion of the individual galaxies in the cluster, and evolve 
the overall system down to the present day, using the LF 
as a constraint a posteriori. This approach also allows, in 
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principle, a more self-consistent description of the chemi- 
cal evolution of the cluster and of its galactic population 
as a whole. 

An improved modelling of the evolution of the cluster, 
taking into account that its galaxies may form at differ- 
ent redshifts, has been introduced bv lChiosil fcOOObh . who 
replaced the usu al integration over the L F with an inte- 
gration over the IPress fc Schechteil (jl974h mass function 
at different redshifts. On the same line, we developed a 
global, self-consistent chemical model for the cluster as a 
whole, following the simultaneous evolution of all of its 
components: the galaxies, the primordial gas, and the gas 
processed and re-ejected via GWs. The approach is also 
somewhat similar in spirit to the "cosmic chemical evo- 
lution" model bv lPei fc Fall (1995j), aimed at calculating 
the global evolution of large volumes of the Universe, pop- 
ulated by a variety of star forming objects. 

Our chemical model for clusters is developed in anal- 
ogy with the usual chemical models for gala xies. These lat- 
ter are sche matically conceived as follows l|Tinslevlll980l 
|Pagellll997l) : 



1. the gas present in the system (usually starting from 
primordial composition) transforms into stars accord- 
ing to some prescribed Star Formation Rate (SFR); 

2. stars are thus formed, distributed according to the 
adopted IMF; 

3. stars return part of their mass in the form of chem- 
ically enriched gas, according to the so-called stellar 
yields (prescriptions derived from stellar evolution and 
nucleosynthesis) ; 



4. this chemically enriched gas mixes with the surround- 
ing gas, causing the chemical evolution of the overall 
interstellar medium (ISM). 

In a cluster, we are interested to model the chemical evo- 
lution of the ICM, and the objects responsible for its en- 
richment are the galaxies, via GWs. In analogy with the 
above scheme, therefore, our chemical model for the clus- 
ter is conceived as follows: 

1. the primordial gas in the ICM gets consumed in time 
by galaxy formation according to some prescribed 
Galactic Formation Rate (GFR); 

2. at each time (redshift) galaxies form distributed in 
mass according to a Galactic Initial Mass Function 
(GIMF), derived from the Press-Schcchter mass func- 
tion suited to that redshift; 

3. galaxies restitute a fraction of their initial mass in the 
form of chemically enriched GWs, according to the 
adopted galactic models; 

4. this enriched gas mixes with and causes the chemi- 
cal evolution of the overall ICM, which includes the 
amount of primordial gas not yet consumed by galaxy 
formation and the gas re-ejected by galaxies in the 
GWs up to the present age. 

The analogy between the two types of models is displayed 
in the scheme above. 

Model equations parallel those of galactic chemical 
models, with the substitutions SFR -> GFR, IMF -> 
Press-Schechter GIMF, stellar yields -> GW yields. The 
main assumptions at the base of the model are as follows. 

— The model is one-zone, namely the cluster is treated as 
a single uniform compound of gas and galaxies, where 
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the metal abundance of the gas evolves in time but is 
otherwise instantly mixed and homogeneous. 
It is worth commenting here on how sensible one-zone 
models are for structures, like clusters, that in current 
cosmological theories form out of hierarchical accretion 
of subunits. We remind here that a one-zone chemical 
model contains no information about the spatial distri- 
bution of its components. Therefore, at high redshifts 
our "model cluster" can be considered simply as the 
sum of the subunits that will later merge and form it, 
irrespectively of whether the cluster has in fact formed 
or not, as a single bound gravitational structure. The 
chemical model at high redshifts just describes the av- 
erage properties of the sum of the "parent subunits" 
of the cluster. 

— The model is calculated assuming the Instantaneous 
Recycling Approximation (IRA), that is assuming that 
galaxies eject the corresponding GWs instantly, as 
soon as they are formed; this is a reasonable approxi- 
mation as the timescales for the onset of the GW are 
generally short, less than 1 Gyr. The IRA could affect 
predictions at high redshifts, where a time-span of a 
few 10 8 yr corresponds to a sizeable gap in redshift; 
but up to redshift z ~ 1, where observational data on 
ICM abundances are available, the effect is minor. 
Later on, the model might be improved in this respect, 
by taking into account the actual delay between the 
formation of a galaxy and the time when its GW is 
expelled; but as it is always better to start with the 
simplest possible assumptions, we adopt the IRA for 
the time being. 

Notice however that the individual galactic models and 
their GW yields are not calculated in the IRA, but 
with a detailed chemical network taking into account 
the different, finite stellar lifetimes at varying stellar 
mass ( Holl and references therein). 

— iMorettil ( 2001 ) investigated a set of pre-chosen, arbi- 
trary GFR laws, providing some indication as to what 
type of galaxy formation history is suited to reproduce 
the observed LF. An early burst of (dwarf) galaxy for- 
mation is needed to reproduce the steep faint end of 
the LF (see Fig. 01; this burst should be followed by a 
more gentle GFR increasing in time (i.e. at decreasing 
redshift) up to a certain epoch, dropping after that. 
Fig. |H1 shows such an example of suitable galaxy for- 
mation history (GFH) and corresponding LF. 

In the present model, we try to include this kind of be- 
haviour self-consistently, by means of a "double infall" 
scheme. 

— Further simplifying hypotheses regard the history and 
the morphology of the galaxies contributing to the 
chemical evolution of the cluster: in facts we have con- 
sidered only elliptical non-interacting galaxies, since 
just early type galaxies seem to be responsible for th e 
metal enrichment of the ICM fsee lArnaud et afll 992). 
Recent arguments about the baryon fraction in galax- 
ies suggest however that early winds migh t be a gen- 
eral requirement for all objects so in this 
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Fig. 8. An example of Galactic Formation History 
[Gyr -1 ] as a function of redshift and corresponding 
Luminosity Function (observational data by Trentham 
1998). Model parameters are listed in Table [21 

sense the distinction between elliptical or other galax- 
ies might be of minor significance. We do not explic- 
itly take into account the possibility of mergers nor of 
episodic star formation after the GW. 

— As to the cosmological parameters, which enter our 
model mostly through the time redshift relation, 
we adopt a flat ACDM cosmology with fiM=0.3, 
Oa=0.7 and Hq = 65 km sec -1 Mpc -1 , corresponding 
to a present age of the Universe, or Hubble time, 
t H = 14.5 Gyr. 

For the Press-Schechter Galactic IMF (see S15.2J1 we 
usually assume n = 1.5 as the index of the power spec- 
trum. 

5.1. Basic equations 

At any time the cluster is assumed to consist of the fol- 
lowing components: 

— The gas with primordial chemical composition and to- 
tal mass Mg t p(t), out of which galaxies are formed. 

— The gas processed through galaxy formation and re- 
expelled by galaxies at the stage of galactic wind. Its 
total mass is M 9i w(t). The wind-processed gas con- 
tains several chemical species (carbon, oxygen, iron, 
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etc.). We will explicitely consider here only oxygen and 
iron, with total masses denoted by M 9} o and M g ^ e , 
respectively. 

We remark here that "wind-processed gas" in our 
models indicates gas that has been processed through 
galaxies, not necessarily all of it through the stars in 
the galaxies. In fact, the GW is made both of metal 
enriched stellar ejecta and of pristine gas "energized" 
enough by stellar feed-back to leave the galaxy as 
wind, without having ever undergone actual star for- 
mation and stellar nucleosynthesis. This pristine gas in 
the GWs should not be confused with what we label 
as "primordial gas" M S) p in our cluster model, meant 
as gas that has never been involved in the process of 
galaxy formation. 

— The galaxies, whose individual mass is the mass in 
stars and remnants left over after the GW stage - 
which occurs instantaneously as soon as the galaxy 
is formed, in the IRA. The total mass in galaxies is 
denoted by M G (t). 

— The total baryonic mass given by 

M b (t) = M g , P {t) + M BtW (t) + M G {t) 

The Dark Matter component does not enter model equa- 
tions as it does not intervene in the chemical evolution of 
the cluster. 

As anticipated, for the primordial gas whence galax- 
ies form we adopt a (double) infall scheme; this scheme is 
meant to describe not quite the accretion of baryons onto 
the cluster region, but the accretion onto the individual 
galaxies or, more precisely, the rate at which baryons be- 
come available (thr ough cooling etc.) f or galaxy formation. 
As pointed out bv lPei fc Fall l| 19951) . infall (or outflow) 
terms can be used even when dealing with large closed 
volumes, or with the whole Universe, as long as they rep- 
resent gas exchange between the individual star forming 
objects and the surrounding medium. Owing to our dou- 
ble infall hypothesis, M b increases in time by accretion of 
primordial gas following the equation 



dM b 
dt 



dM, 



dt 



= Ate +B(t-t )e ^ (3) 



in J 



The first term represents a first, fast (~ 0.5 Gyr) "infall 
episode" forming dwarf galaxies at high redshift; these ob- 
jects are presumably responsible for reionization and/or 
reheating of the intergalactic medium, so that galaxy for- 
mation is temporarily halted. Later on, gas cools down 
again and galaxy formation progressively sets in again in 
a second, smoother and prolongued "infall phase" repre- 
sented by the second term (active for t > to). This double 
infall scheme was chosen so as to yield a shape of the GFH 
suitable to repr oduce the obse rved LF (from a preliminary 
investigation bv lMorettfeOOll see Fig.|HJ), as well as being 
reminiscent of current evidence and theories about reion- 
ization at high redshift and of the observed cosmic SFH 
of the Universe at lower redshift (the "Madau-plot" ) . 

The double infall scheme causes the baryonic mass of 
the cluster to increase up to a final value M^t at the 



present time iff (Hubble time, corresponding to z = 0). 
We adopt this final total baryonic mass M b ,T as the nor- 
malization mass to which all values are scaled: 

M b (t H ) = M b>T = 1 (4) 

Hence, in Eq. J3J we can tune the timescales t% and Ti of 
the two infall episodes, as well as the amount of mass in- 
volved in the second infall phase (set by the parameter B), 
while the parameter A is fixed by the normalization 

t H~ t a ' 

(T2 + iff - t ) e T 2 



1-Bt 2 



A = 



n 



n 



(n + t H ) e T1 



With the adopted normalization the masses of the 
different cluster components (M 9j p, M 9t w and M G ) effec- 
tively become, as well as M b , mass fractions of the final 
baryonic mass M b: r, and as such they are dimensionless 
quantities. The equations governing their temporal evolu- 
tion are: 



dM, 



9,P 



dt 



-*g(*) 



dM, 



g,P 



dt 



dM, 



dt 



dMc 
dt 



* G (t) [1 - Y G (t)} 



(5) 



(6) 



(7) 



The equations governing the masses of oxygen and iron 
M 9j o and M 9j p e expelled into the ICM by GWs are: 

dA/ fl ,o 



di 



dM, 



9,Fe 



dt 



^ G {t)Y G ,o{t) 



*G(t)Y GjFe (t) 



(8) 



(9) 



In the equations above, ^> G represents the Galaxy 
Formation Rate (GFR) — the analogue of the SFR in 
classical chemical models of the interstellar gas — while 
Yq , Y G ,o and Y G ,Fe are the "yields" per galactic genera- 
tion of total gas, oxygen, and iron respectively; all these 
quantities will be specified in the following sections. 
The total gas mass in the ICM is: 

M g = Mg, P + M,,u 

and its oxygen and iron abundances are given by: 



X = 



M, 



g,0 



M n 



Xfc — 



M, 



9,Fe 



M n 



lg IVXg 

Finally, the initial conditions for our set of equations are 

M s> p(0) = M giW (0)=0 M G (0)=0 (10) 

together with 

M 9 , o (0) = M g , Fe {0) = (11) 

The set of equations is integrated with a 4 t7l -order Runge- 
Kutta method, dividing the integration time-span [0 — iff] 
into 1000 timesteps. The consistency of the solution is 
secured by the fact that at the final epoch i = iff the sum 
of the different cluster components is: 

M fl (iff) + M G (iff) = M 6 , T = l 
within a small error (of the order of 10 -3 ). 
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5.2. Galactic Initial Mass Function 

The GIMF adopted here stems from the Press-Schechter 
(PS) mass function, originally devised to describe the 
collapse and formation of dark matter haloes from the 
primordial spectrum of perturbations. Following e.g. 
lLacev fc Cote dl994|) . the mass function of collapsed ob- 
jects at redshift z is: 



dhiM 



M 



x exp 



M 



M*(z) 



(12) 



where n is the index of the cosmological power spectrum, 
z is the redshift corresponding to epoch t, and M* is the 
cut-off mass, or characteristic mass, at that redshift: 



M*(z) = M Nor x {1 + z)- — , M Nor = M*{z 



0) 
(13) 

Eq. I|12|) yields the fraction of mass in dark matter haloes 
with mass M at redshift z{t), per logarithmic mass inter- 
val. As discussed at the beginning of for the purpose of 
the present study we assume that the same formula pro- 
vides also the mass spectrum of the galaxies forming at 
redshift z{t). For a fixed baryonic fraction in the Universe, 
typically f b = 0.1, the ratio M/M*(z) in Eq. ijT2"|> can be 
considered both as the ratio between a dark halo mass 
and the corresponding "characteristic mass" at that red- 
shift, and as the ratio between the corresponding bary- 
onic masses. Namely, as the PS mass function is expressed 
solely in term of a mass ratio, it can be used both as a 
mass function of the dark haloes, and as a mass function of 
the corresponding baryonic components. In the latter case, 
M*(z) becomes the baryonic mass of the (proto)galaxy 
hosted by the "characteristic halo" at redshift z, or the 
typical mass-scale of (proto)galaxies at that redshift. In 
this context, M]\f or is the mass of baryons contained in 
the typical dark halo mass collapsing at the present-day. 

In our model the GIMF <fi G is defined as the number 
of galaxies per mass interval, given by: 



dfW)} 



M 2 d\nM 



(14) 



where the time dependence is due to M*(z(t)) on the 
right-hand side. 

While Eqs. (H2J and (O define the shape of the GIMF, 
we need to define also the mass range [M i: M u ] within 
which galaxies can form at any redshift; in principle these 
mass limits can vary with redshift. 

For the upper and lower mass limits we set, respec- 
tively: 



M u (z) = 7 M*(z) 



M u {z) 



M.(z) 
7 



with 7 ~ 100; that is to say, at any given redshift galax- 
ies form with masses within two orders of magnitude the 
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Fig. 9. Evolution of the characteristic mass M* (in this 
example, with Mnov = 3x 1O 13 M0) and of the minimum 
and maximum mass limits for the GIMF, Mi = M*/100 
and M u = 100 M*. 

characteristic mass at that redshift. Fig. [U| shows an ex- 
ample of the evolution of M* , Mi and M u as a function of 
redshift (example with M Nor = 3 x 10 13 M and 7 = 100). 

The resulting galactic mass function is finally normal- 
ized in mass, at any redshift, according to: 



M„(t) 



M(j) G {M,t) dM = 1 



M z (t) 



so that <pc(M,t) properly becomes the distribution func- 
tion of the masses of the galaxies formed at epoch t. 

5.3. Global galactic yields 

Once the mass function of galaxies and its mass limits 
are specified as a function of time (redshift), the global 
galactic yield of gas and metals of each galactic generation 
can be determined out of integration of the GW yields of 
individual galaxies over the mass function. 

Our galactic models, discussed in W6.2\ provide the 
mass of gas Eqw(M, t), iron Ecw,Fe(M, t) and oxygen 
Egw,o(M 7 t), ejected in the GW by a galaxy of given ini- 
tial mass M and formation epoch t. The global quanti- 
ties of gas, iron and oxygen instantaneously recycled and 
re-ejected into the ICM by an entire galaxy generation 
formed at epoch t can be determined as: 



Y G (t)= / E GW (M, t) G (M, t) dM 

JMi(t) 



(15) 



rM u (t) 

Y G>Fe (t)= E GWiFe (M,t) 4> G (M,t)dM (16) 



/■M„(t) 

Y G .o(t)= E GWt0 (M,t) d> G (M,t)dM (17) 

JMi(t) 
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These are the quantities that enter Eqs. i|til9|) and which 
must be evaluated at each timestep t of the model, since all 
the quantities involved (GIMF, mass limits, GW ejecta) 
depend on the epoch t. The mass that remains locked in 
the galaxies ("galactic remnants", see S13.2J) per galactic 
generation is: 

R G (t)= / R gal (M,t)<p G (M,t)dM = l-Y G (t) 

JMi(i) 

(18) 

entering Eq. 10. 

5.4. Galactic Formation Rate 

As to the Gal actic Formation Rate, we adopt a simple 
ISchmidtl ||1 9501) law 

^ G = uM g . P (19) 

where galaxies are assumed to form solely from the pri- 
mordial gas component with a rate that is directly pro- 
portional to the available "fuel" . With this law, the time 
evolution of the GFR is primarily set by the "infall his- 
tory" of the primordial gas. 

The Galaxy Formation History is supposed to start at 
epoch z = 20 (^ G (t) = for z > 20), and it is halted when 
the characteristic mass forming at that redshift, M*(z), 
becomes larger than that of the brightest observed galaxy 
in the LF. In our simple "hierarchical approach" , namely, 
the formation of the brightest /largest galaxy corresponds 
to the lowest typical redshift at which galaxy formation is 
active. 

6. A fiducial model 

The main observational constraint to calibrate the Galaxy 
Formation History is t he ob served luminosity function. 
The LF by iTrenthaml l(l998L Fig. 0) shows the charac- 
teristic Schechter-like exponential cut-off at bright mag- 
nitudes (Mb < —20), but it also steepens at the faint end 
(Mb > —14). This large number of faint galaxies seems to 
require a prominent episode of for mation of small galax- 
ies at high redshift l|Morettil 12001 and Fig. 0, a feature 
which inspired our present "double infall" prescription. 

The LF of cluster galaxies in our models is calculated 
by assigning to each galaxy of "remnant" R ga i(M, Zf or ) 
the relevant M/L ratio for the age corresponding to its 
redshift of formation Zf or (see the discussion in H3.3|l . 

As to the normalization of the total number of galaxies 
in the cluster, our models are calibrated also so that the 
final proportion of intra-cluster gas and galaxies respects 
the observational constraint of the ICMLR^37 f iJ4.2|) . 
within a 20%. 

Here we present and discuss in detail our "fiducial 
model", Model A whose parameters and final results are 
listed in Table PJ1 In this model we adopt as the character- 
istic mass for the GIMF ( $57% 

M Nor = M*(z = 0) = 3 x 1O 13 M 




Fig. 10. Luminosity Function of cluster galaxies as pre- 
dicted by our "fid ucial" Model A; triangles are the ob- 
servational data bv ITrenthaml 1|1998(1 . Both the model LF 
and data are shown normalized to the same total (unit) 
number of objects in the observed range (Mb between -22 
and -11 mag). 
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Fig. 11. Top panel: GFH of our Model A as a function of 
redshift; *5f G is expressed in mass fraction of the cluster per 
Gyr. Bottom panel: comparison of th e shape of ou r GFH 
with the "Madau plot" (data from ISteidel et alJ ll999): 
the theoretical curve ^Sq is arbitrarily normalized so that 
its maximum value reaches log( v I'G) =_ 0.8, similar to the 
maximum in the data. 

which is of the orde r of the characteristi c mass of galaxy 
clusters nowadays (jGirardi et alJ Il998j) once a factor 
fb = 0.1 is applied to scale the total mass to the mass 
of the sole baryonic component. 

With the adopted parameters the final ICMLR~33, 
the observed LF is reasonably well mimicked (Fig. I10|) 
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0.09 
0.35 


1.50 
1.25 


0.0184 
0.0215 


8.2 
9.4 


2.26 
2.29 


40.6 
39.2 


0.247 
0.299 


0.163 
0.164 


Fig.H 


PNJ 




3.el3 


input GFH 


as in 


Fig-El 


0.0061 


5.8 


1.56 


30.7 


0.108 


-0.060 



Table 2. Parameters and results of the cluster models 
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Fig. 12. Time evolution of the cluster components for 
Model A 

and the galactic GFH proceeds until z ~1 (Fig. 111(1 : 
at lower redshifts, the typical galaxies forming would be 
more massive/bright than observed in the LF. The steep, 
faint end of the LF is reproduced by means of an early 
galaxy formation activity connected to the first infall 
episode (Fig. 1111 top panel); this early activity is short 
(ri ~ 0.15 Gyr) and can be identified with the forma- 
tion of the first objects responsible for reionization. In 
fact, the second phase of galaxy formation (correspond- 
ing to the "second infall") slowly sets in starting from 
to = 0.8 Gyr, in broad agreement with recent observa- 
tions suggesti ng that the reionization/reheating era too k 
place at z > 6 ijDiorgovski et al.l200li [Becker et al.l200l|) . 
Hence, we can loosely associate the trough in the galaxy 
formation activity at the transition between the two "infall 
episodes" (Fig. ITT)) with the reheating of the intergalactic 
medium following the formation of the first galactic ob- 



■ Mushotzky & Loewenstein (1997) 
• Matsumoto et al. (2000) 
A Fukazawa et al. (1998) 




1 2 3456789 10 

(1+z) 

Fig. 13. Evolution of the iron abundance in the ICM. At 
high rcdshift [z > 1.5) metallicities for Damped Lyman a 
systems and Lyman Break Galaxies are shown for the sake 
of qualitative comparison ( metallicity ra nges for LBGs 
and DLAs at z = 3 are from |Pettinil l20001. 

jects; the subsequent slow cooling down of the gas fuels 
the GFR at z < 6. The GFH eventually halts so abruptly 
due to the simplified prescriptions in the model ( U5.4jl : we 
remind also that only early-type galaxies (E/S0) are con- 
sidered in our cluster model, while in the Madau plot for 
field galaxies, most of the SF at z < 1 is due to spirals. 

In Fig. EH the dotted part of the histogram represent 
the predicted extension of the LF to the range of galax- 
ies fainter than the observational limit. Although these 
small, faint galaxies dominate in number, their contribu- 
tion in terms of luminosity o r stellar mass is a negligible 
fraction of the total (see also lThomasl ll999 1 : in the Local 
Group as well, the numerous dwarf galaxies play a neg- 
ligible role in the global mass and luminosity budget. In 
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clusters, these faint objects might also have been disrupted 
and be nowada ys dispersed as a diffuse intra-cluster stel- 
lar component (iGregg fc Westlll998tlFerguson et al,lll99^: 
lArnaboldi et alJl2002i and references therein). 

Fig. El shows the evolution of the mass fractions in 
the various cluster components: the mass in primordial 
gas, in wind-processed gas and in galaxies. The "double 
infall" pattern can be recognized in the evolution of the 
primordial gas component; within the first Gyr, the early 
infall phase produces a minor peak — first the gas mass 
increases due to infall, then it tends to decrease, being 
consumed by galaxy formation; later on, again the increase 
due to infall competes with the consumption by galaxy 
formation, until the latter halts at z — 1 . At the end of the 
evolution, wind-processed gas is ~1.5 times the mass in 
galaxies, in agreement with qualitative expectations from 
individual galaxies with the PNJ IMF (323- The total 
mass in gas is about 6 times the mass in galaxies; about 
70% of the ICM is primordial gas never involved in the 
galaxy formation process. 

Fig. El shows the time evolution of metallicity (iron 
abundance) in the ICM; it peaks at z ~ 1, when the most 
massive galaxies are formed and release their metal pro- 
duction, then it decreases because primordial gas is still 
being accreted onto the cluster, following our infall pre- 
scription, even after galaxy formation and metal produc- 
tion are over. The final metallicity is low by a factor of 3 
with respect to observations, in agreement with the low 
overall IMLR (~0.006, Table 0). This result is easily un- 
derstood since most galaxies (in terms of mass involved) 
form at redshifts z < 5, with a low corresponding IMLR 
(see the discussion in 21 and Table |2J). 

The predicted [O/Fe] ratio is roughly solar (Table 0). 
If we were to compensate for the low iron production by 
increasing the rate and iron contribution of SN la by the 
required factor of 3, the resulting [O/Fe] ratio would de- 
crease by 0.5 dex, and become too low with respect to 
observations, which indicate rati os between solar and su - 
persolar for the bulk of the ICM IjFinoguenov et all 2000). 

6.1. Comparison to the Salpeter case 

For the sake of comparison, we calculate an analogous 
model (Model AS, with parameters listed in Table 01 
adopting GW ejecta and "galactic remnants" as computed 
with the Salpeter IMF. Fig. [His the analogue of Figs.ITT)! 
to El for Model AS. The LF is reproduced with a GFH 
qualitatively similar to that of Model A, but both the 
metal abundances in the ICM and the IMLR (Table 
are far from observational data. Besides, in this case the 
mass of wind-ejected gas is comparable to, or lower (70%) 
than, the mass in galaxies and constitutes a minor fraction 
(13%) of the total final ICM. 

Namely, with a Salpeter IMF galaxies produce by far 
too little metals per mass stored in galaxies. With the PNJ 
IMF things visibly improve (Fig. I13J) . still the expected 
metallicity is ~3 times lower than observed. 



7. Varying the Galaxy Formation History 

In this Section we explore alternative models with 
"slower" or "faster" GFHs, by considering different char- 
acteristic masses M^ or for the GIMF fi l5.2[) . Qualitatively, 
increasing M^ or , the galaxies forming at a given redshift 
are intrinsically more massive and luminous, so that to 
reproduce the observed LF the optimal GFH is "antici- 
pated" , or skewed to higher z. Thus, M^ or is linked to the 
peak (and halt) of the corresponding GFH. Similar effects 
could be obtained also by changing the spectral index n, 
for a shallower power-spectrum also induces the formation 
of intrinsically more massive galaxies at a given redshift, 
henceforth "anticipating" the required GFH. We explored 
models with n values different from the -1.5 adopted here, 
but we do not report them here for the sake of brevity, as 
conclusions would be similar to those with varying M^ or . 

Models Ml, M2 and M3 are analogous to Model A, 
but for different values of M^ or ; in each model, the "sec- 
ondary parameters" governing the GFH (£?, t\, T2, v) are 
optimized so as to reproduce the observed LF and ICMLR. 
Model Ml has M Nor = 5 x 10 12 M Q , comparable to the 
characteristic mass at z = in the field, which in current 
cosmologies is of the order of a few 10 13 M Q (total mass, 
dark matter included). Models M2 and M3 have a value 
of Mjvor larger than in Model A; their Mjvor corresponds 
to the typical baryonic content of rich clusters, with total 
mass of the order of 10 15 M@. 

As anticipated, at increasing M^or (from Ml to M3), 
the GFH is skewed to higher z (Fig. El d-panels); also 
the bulk of the metal production is correspondingly an- 
ticipated. The main consequence is that the global IMLR 
increase's with M Nor (TablcEJ, since a given L b bin in the 
LF now corresponds to older, dimmer and more massive 
galaxies, which eject both more gas and more metals into 
the ICM, and have a higher characteristic IMLRgvk ^3.21 
and 13.3(1 . Nevertheless, even in model M3 corresponding 
to the largest Mjvor, the predicted IMLR is still low by 
a factor of 2 with respect to observations. This is due to 
the fact that in all models most galaxies in terms of mass 
involved form anyways at z <5, with a low characteristic 
IMLR, as was the case for model A. Correspondingly, the 
predicted metallicities are too low. 

Hence, considering different GFHs — in particular, 
"anticipated" GFHs with respect to model A, as in mod- 
els M2 and M3 — improves the situation with metal pro- 
duction, but only slightly. The conclusion in the previous 
section still holds, that if an IMF varying with redshift 
performs much better than the Salpeter models, the metal 
production provided by the PNJ models is still low. 

8. Choosing the IMF "ad hoc" 

In the previous Sections we showed that a varying IMF like 
the PNJ IMF, skewed toward more massive stars (or more 
precisely, having a lower locked-up fraction) at higher red- 
shifts, yields much better results for the predicted ICM 
metallicity, than the standard Salpeter IMF. Nevertheless, 
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Fig. 14. Same as Figs. ITUHT31 but for Model AS based on the Salpeter IMF. 



a still larger metal production seems to be required, by a 
factor of 2-3, than obtained with our PNJ models. In these 
models, aside from the variation with galaxy mass due to 
density effects, the peak mass of the IMF mainly varies 
with rcdshift because of the CMB temperature f ^3.1[) : ap- 
parently, a larger variation is required to reach the re- 
quired metal production. I nterestingly, this is in line wit h 
other recent results: both iHernandez fc Ferraral l|200ll) . 
from the numb er of metal-poor stars i n the halo of the 
Milky Way, and lFinoguenov et al.1 l|2003h . from arguments 
related to the different yield in groups and clusters, sug- 
gest a typical peak-mass of the IMF increasing with red- 
shift faster than what expected from the pure Jeans effect 
of the CMB temperature. 

We consider in this section some cluster models as- 
suming an IMF that is more efficient in terms of gas and 
metal ejection in the galactic wind. We achieve this by 
"fixing" the galaxy models to be the PNJ models with 
zj or =15: the GFH in the cluster still extends in time down 
to z ~1 (or lower), but now we assume that galaxies at 
any redshifts form with the properties (galactic yields and 
remnants) of the PNJ models with z/ or =15. Namely, the 
integrals IT51 through ITS1 in tJ5.3l are now computed always 
from the models with Zf or —15, rather than with a Zf or 



running with the actual epoch z of cluster evolution. This 
is a simple artifact to mantain an IMF "more efficient" 
than the CMB-regulated PNJ IMF. Luminosities are of 
course computed for the actual age t(z) of the galaxies. 
The choice of the models with z/ or =15 is ad hoc to ob- 
tain the correct final IMLR and metallicity (see 21 an d 
Fig.UH). The IMLR resulting from the "hierarchical" GFH 
is somewhat lower than the monolithic case with the same 
z for — 15 because, with galaxy formation extending down 
to redshift 1-2, galaxies are on average younger and more 
luminous than in the monolithic case, lowering the typical 
IMLRgvk f° r a given mass of ejected iron. 

In Fig. El w e show the cluster models AfY and M2fY, 
analogous to models A and M2 (i.e. with the same M^or = 
3 x 10 13 and 10 14 , respectively) but calculated with "fixed 
galactic yields" as described above; parameters and results 
are listed in Table [21 In both cases the final metallicity 
and IMLR are in good agreement with observations, the 
mass in ejected wind is ~2.3 times the mass in galaxies 
and constitutes roughly half of the whole ICM gas. Notice 
also that, for these models reproducing the correct final 
ICM metallicity, the predicted metallicity at redshift 3- 
5 falls in the range of Lyman Break Galaxies, which are 
in fact considered to trace the high redshift counterpart 
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Fig. 15. Results for Models Ml, M2 and M3 which differ from our fiducial Model A in the characteristic mass Mm or, 
increasing from Ml to M3 (Tabled- (a) Luminosity Function; (b) evolution of the iron abundance; (c) evolution of 
the cluster components; (d) Galactic Formation History. Lines, symbols and data as in Figs. I1UI to 1131 



of present-day massive galaxies in high density regions of 
the Universe, such as clusters. 

In model AfY, the GFH peaks at z ~1 (similarly to 
the Madau plot), and correspondingly a noticeable peak 
in metallicity is predicted at the same redshift, for no fur- 
ther enrichment takes place for the primordial gas infalling 



onto the cluster later on. This fast metallicity evolution is 
not really observed — although very few data points are 
available for redshifts z > 0.6, where most of the evolution 
takes place. 

A milder behaviour, somewhat more compatible with 
data, is found for model M2fY, having a characteristic 
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Fig. 16. Results for Models AfYzl5 and M2fYzl5, analogous to models A and M2 but with fixed galactic yields 
corresponding to the PNJ models with z = 15. Panels and symbols as in Fig. 1151 



z = mass Mjvor of the order of the baryonic mass in rich 
clusters. For this model, the GFH peaks at z ~2 and halts 
at z ~1.5, the metallicity peaks at z ~1.5 while evolving 
rather smoothly out to z ~1. We consider model M2fY as 
our best model. 



9. Abundance ratios in the ICM 

As mentioned in the introduction, the typical [a/Fe] abun- 
dance ratios in the ICM are not firmly established from 
observations. With the first ASCA data, it was debated 
if the ICM is typically a-enrich ed or if its abundances 
are consistent with solar ratios l|Mushotzkv et all Il996t 
llshimaru fc Arimotoi ri997). Recent studies have revealed 
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Fig. 17. Evolution of the [O/Fe] abundance ratio as a 
function of redshift for our Models A, AS and AfY, com- 
pared to observations for low redshift clusters 



that abundance ratios may change in different regions of 
the cluster, from roughly solar in the central regions t o 
a-enriched in the outer regions l|Finoguenov et al]l2000l) . 
This might suggest competing, or rather overlapping, ex- 
traction mechanisms of metals from galaxies: SN II ejecta 
would dominate in GWs, while additional SN la prod- 
ucts would be extracted from the galaxies due to ram- 
pressure stripping, mor e effective in the central regions 



( Finoguc nov et al 



g, more 
1. 2000). 



Our model at present does not deal with gradients of 
abundance or of abundance ratios in the ICM, however 
it provides the global oxygen mass produced by galaxies 
and ejected into the ICM by GWs, so that we can estimate 
the typical average [O/Fe] ratio in the ICM, with oxygen 
being the best tracer of a-elements. In Fig.Elwe show the 
evolution of the [O /Fe] ratio in the ICM as predicted from 
our reference Model A; [O/Fe] predictions for models Ml- 
M2-M3 are very similar so these models are not shown in 
the plot. [O/Fe] is supersolar at very high redshifts (where 
the PNJ IMF favours more massive stars and hence SN II) , 
decreasing down to marginally subsolar values. Model AS 
with the Salpeter IMF is also shown for comparison; the 
final [O/Fe] value in very close to that of Model A. 

Also shown is model AfY with the "fixed yields" PNJ 
models (model M2fY is basically indistinguishable from 
AfY in this plot and is not shown); the predicted final 
[O/Fe] is slightly supersolar (around +0.2 dex). Due to the 
"fixed yields" assumption in this model, [O /Fe] is roughly 
constant throughout the evolution and its value is charac- 
teristic of the PNJ models with Zf or — 15 (see 2J. 

Given all the above mentioned caveats about the com- 
plexity and uncertainties of the empirical evidence, we 
consider all these final values to be in broad agreement 
with observations. Notice however that all models predict 
a very slow evolution (if any) in the [O/Fe] ratio, so that 
no sistematic trend is expected to be seen in for this ob- 
servable until z > 3 at the earliest. 



10. Summary and conclusions 

Galactic winds from elliptical galaxies are the most likely 
source of the chemical enrichment in the ICM. In this 
scenario, various studies in literature suggest that a non- 
standard IMF must be invoked for elliptical galaxies, if 
we are to a c count for the metallic ity of the ICM (Sj2J. 
IChiosi et all ljl99Sl) . IChios] l|2000bh calculated models of 
elliptical galaxies adopting the variable IMF by PNJ, 
whose behaviour in time and space is sensitive to the 
physical conditions of the star-forming gas; in particular, 
the low-mass cut behaves as a sort of Jeans mass (JU- 
This IMF naturally predicts a lower locked-up fraction in 
the early galactic stages, especially for massive ellipticals 
and/or for high redshifts of formation — the latter feature 
being related to the increasing temperature of the CMB. 
These galactic models successfully re produce a variety o f 
observational properties of ellipticals l|Chiosi et al.l fl998') . 
In H3.2l we show how models calculated with the PNJ IMF 
predict galaxies to eject both more metals and more gas 
than "standard" models based on the Salpeter IMF. We 
also show that, for the more massive objects that play the 
main role in the mass and metal enrichment budget of the 
cluster, the characteristic IMLR of individual galaxies is 
much higher for models with the PNJ IMF than with the 
Salpeter IMF (Q. 

To assess the effect of these new galactic models on 
the ICM, first we assumed all galaxies to be coeval and 
we integrated the corresponding ejecta over the observed 
luminosity function, as in the standard "monolithic ap- 
proach" (2J. The Salpeter models fail in reproducing the 
observed IMLR in clusters by an order of magnitude, while 
better results are obtained with the the PNJ IMF, espe- 
cially with the models corresponding to a redshifts of for- 
mation Zf or ~13. We also discuss the evidence of dilution 
of the galactic wind ejecta with primordial gas, introduc- 
ing the concept of the intra-cluster-mass-to-light ratio 
(ICMLR); for the typical luminosities and M/L ratios of 
our model galaxies, a dilution of a factor 2-3 for the PNJ 
models (a factor of 6 for the Salpeter models) is necessary 
to recover the typical amount of gas vs. galaxies observed 
in rich clusters. 

Alternatively to the standard monolithic approach, we 
considered a simple hierarchical picture in which galaxy 
formation extends in time, with a characteristic galaxy 
mass increasing at decreasing redshift, as from Press- 
Schechter theory. We developed a toy-model following 
the chemical evolution of the ICM in connection with the 
galaxy formation history of cluster galaxies (fJSJ). The GFH 
is calibrated so as to reproduce the observed present-day 
galactic LF in clusters, and the observed amount of ICM 
gas (in terms of ICMLR) . This seems to require two phases 
of galaxy formation: an early phase forming dwarf galax- 
ies populating the steep, faint tail of the LF, followed by 
a more gentle GF activity peaking at z ~ 1 — 2. This sce- 
nario is in broad agreement with (a) present theories and 
observational evidence about reionization and reheating 
of the Universe at z > 6, and (b) the observed trend of 
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cosmic star formation history — although in our cluster 
models the GFH is "anticipated" and halts earlier with 
respect to the field-based Madau-plot, since we are con- 
sidering early-type galaxies. We mimick this "bimodal" 
galaxy formation history by means of a double infall pre- 
scription in our chemical model for clusters. 

A satisfactory match with the observed LF can be 
obtained both using galactic models with the PNJ IMF 
and with the Salpeter IMF. However, models with the 
PNJ IMF, besides bei ng favoured on th e base of their 
photometric properties ijChiosi et al.lll998j) . provide much 
improved predictions about the metal enrichment of the 
ICM; Salpeter-based models fail in this respect (SjlJ). Still, 
the PNJ models are short of metal production by a factor 
of 2-3, and the conclusion does not depend much on the 
details of the GFH (iQ. In fact, with such extended GFHs 
the bulk of galaxies (in terms of mass involved) form at 
redshifts z < 5, where the PNJ models do have a typical 
IMLR which is a factor of 2-3 too low (©. 

The discrepancy cannot be cured simply by increas- 
ing the (quite uncertain) rate of SN la in the models, 
since the resulting [O /Fe] ratio would correspondingly de- 
crease by 0.3-0.5 dex at odds with observations. Our mod- 
els presently give the correct relative contribution of SN la 
and SN II as demonstrated by the predicted [O/Fe] ratio 
in the ICM ( solar or slightly supersolar, and strong 
variations of the SN la rate are not allowed within the 
constraint of the observed relative abundances. 

The implication is that an IMF with a more extreme 
behaviour than assumed in our galactic models is needed. 
Variations of the stellar Jeans mass with redshift in- 
duced by the CMB temperature appear to be too mild 
to produce the observed metal c ontent in clusters, in 
agreement with other recent r esults ijHernandez fe Ferraral 
200ll lFinoguenov et al.l2003f) . In fact, from the monolithic 
approach an IMF is required, behaving like that in the 
PNJ models with Zf or ~13-15. 

We thus computed some cluster models with the IMF 
chosen "ad hoc" to induce a metal ejection large enough to 
enrich the ICM to the observed abundance level ( 3HJ - Our 
best model is M2fY (Fig.HBJ), with a GFH peaking at z ~2 
and halted by z ~1.5. This models predicts a peak in the 
metallicity of the ICM around z ~1.5, decreasing at higher 
redshifts; at z=3-4 the typical metallicity corresponds to 
observed Lyman-break galaxies. In this model, the wind- 
ejected gas is more than twice the mass in galaxies, and 
constitutes roughly half of the ICM mass. 

In general, when galactic models with the PNJ IMF 
are adopted, the mass globally ejected in GWs exceeds the 
mass stored in stars by a factor of 1.5-2 — arguments in 
favour of large outflows of baryons from the main body o f 
galaxies have been recently advanced also bv ISilkl l|2002h . 
Although this is not enough to account for the total mass 
in the ICM, it is nevertheless a non-negligible fraction 
of the intra cluster gas, especially for less rich clusters. 
This suggestion is reinforced by recent re sults obtained 
with dynamical models of gal axy formation ijCarraro et alJ 
l200lt IChiosi fc Carrarol 2002) , showing that the process of 



galaxy formation is so unefficient that roughly 75% of the 
gas initially available is "wasted" and expelled, while only 
a 25% remains locked in the final galaxy. These dynamical 
models, even adopting a constant IMF, suggest a typical 
ratio of 1:3 between the matter stored in galaxies and that 
re-ejected in the inter-galactic, or intra-cluster, medium. 
It is suggestive to speculate that 3:1 or so could set a 
minimal value for the gas-to-galaxy mass ratio in clusters; 
this value could then be typical of poor clusters, whose 
ICM would be dominated by gas of galactic origin. This 
is also in line with the observational evidence that the 
energetics of the ICM in the less rich and luminous clusters 
are dominated by non-gravitational effects, most likely 
related to the energy feed-back from cluster galaxies. 
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